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ABSTRACT 

We report the discovery of a one-sided 3.6" (24 kpc, projected) long jet in the high-redshift, z=4.72, 
quasar GB 1428+4217 in new Chandra X-ray and VLA radio observations. This is the highest redshift 
kiloparsec-scale X-ray/radio jet known. Analysis of archival VLBI 2.3 and 8.6 GHz data reveal a 
faint one-sided jet extending out to ^200 parsecs and aligned to within ^30° of the Chandra /VLA 
emission. The 3.6" distant knot is not detected in an archival HST image, and its broad-band spectral 
energy distribution is consistent with an origin from inverse Compton scattering of cosmic microwave 
background photons for the X-rays. Assuming also equipartition between the radiating particles and 
magnetic field, the implied jet Lorentz factor is ~ 5. This is similar to the other two known z ~ 4 kpc- 
scale X-ray jet cases and smaller than typically inferred in lower-redshift cases. Although there are still 
but a few such very high-redshift quasar X-ray jets known, for an inverse Compton origin, the present 
data suggest that they are less relativistic on large-scales than their lower-redshift counterparts. 

Subject headings: Galaxies: active — galaxies: jets — quasars: individual (GB 1428+4217) — radia- 
tion mechanisms: non-thermal — radio continuum: galaxies — X-rays: galaxies 



1. INTRODUCTION 

The z = 4.72 quasar GB 1428+4217 (B3 1428+422) 
was identified in a search for high-redshift objects 
through targete d optical spectroscopy of flat-spectrum 
radio sources (|Hook fe McMahonl 11998] : IFabian et all 
fl997h . It is a luminous X-ray source with detected 
X-ray and radio va r iability charac t eristic of a blazar 
(IFabian et all 119991 : iWorslev et ail 120061: I Veres et al.l 
2010). On parsec-scale, VLBA 15 GHz images indicate 
a dominant high brightness temperature, 71, ~ (4 — 6) x 
10 11 K, core component w ith a faint one-sided jet-like ex- 
tension (|Veres et al.ll2010L and references therein). The 
high luminosity, variability, and radio compactness are 
all properties consistent with Doppler beaming of emis- 
sion from a relativistic jet aligned close to our line of 
sight ( see IFabian et al.lll999D . 

Such high-redshift radio/X-ray sources offer a unique 
glimpse into powerful outbursts from active galactic nu- 
clei (AGN) in the early Universe. As the ambient 
medium into which large-scale jets propagate (including 
host galaxy environments and intergalactic medium) is 
expe cted to be drastically different at such earl y epochs 
(e.g. JDe You^l200l iMilev fc De Breuckll2008D . studies 
of large-scale jet structures allow us to probe radio source 
interactions with their environment. Motivated by X- 
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ray detections of kpc-scale jets in t wo very high-redshift 
quasars, 1745+624 at z = 3.9 ({Cheung et al.l I2006D 
and GB 1508+57 1 4 at z = 4.3 (jSiemiginowska et all 
120031: lYuan et all 120031 ). and the noticeable dearth 
of Chandra observations of z ^ 2 jet systems (e.g. 



Kataoka fc Stawar9 120051 lHarris fc Krawczvnskil 120061 : 
Massaro et alj 120111 ). we began a program to obtain 



arcsecond-res qlution radio and X-ray i maging of more 
such systems (jCheung et alJl20~0 5. 2008) with the aim to 
understand the physics of the highest-redshift relativistic 
jets. 

Using Chandra X-ray Observatory and NRACQ 
Very Large Array (VLA) imaging observations of 
GB 1428+4217, we discovered an X-ray/radio jet sep- 
arated from the nucleus by 3.6" (24 kpc, projected^. At 
z = 4.72, this is the most distant kpc-scale jet imaged 
in X-rays. No significant optical emission is detected 
from the 3.6" knot in an archival Hubble Space Telescope 
(HST) image. To probe smaller scale emission, we im- 
aged archival very long baseline interferometry (VLBI) 
2.3 and 8.6 GHz data, revealing a one-sided ~200 par- 
sec long jet aligned within ~ 30° of the kpc-scale Chan- 
dra /VLA structure. In the following, we present these 
multi- wavelength observations (Section 2), and go on to 
discuss the physical parameters of the large-scale outflow 
in terms of inverse Compton emission models, comparing 
this case to other X-ray detected AGN jets (Section 3). 

2. RADIO, X-RAY, AND OPTICAL OBSERVATIONS 

As part of a larger survey search for kpc-scale jets 
in a sample of z £ 3.4 flat-spectrum radio sources 
(|Cheung et al.ll2005D . we obtained VLA A -array observa- 

7 The National Radio Astronomy Observatory is operated by 
Associated Universities, Inc. under a cooperative agreement with 
the National Science Foundation. 

8 Adopting H = 71 km s" 1 Mpc -1 , Q M = 0.27 and U A = 0.73, 
the quasar is at a luminosity distance, Dl = 44.5 Gpc and 1"= 
6.59 kpc. 



2 



Cheung et al. 




Fig. 1— Chandra 0.3 - 7 keV (color) and VLA 1.4 GHz (con- 
tours) images of GB 1428+4217 showing the bright core and faint 
~ 3.6" distant jet knot at PA ~ 295°. Coordinates are in J2000.0 
equinox. The X-ray data are binned by 1/2 of the native 0.492" 
pixels and Gaussian smoothed with kernel radius of 3 pixels. The 
10 radio contours start at 0.17 mjy/bm (4 times the off-source 
rms) increasing by factors of 2 up to 87 mjy/bm (peak is 155.4 
mjy/bm) with circular beamsize = 1.5". 

tions of GB 1428+4217 on 2004 Dec 6 (program AC755) 
at 1.4 and 4.9 GHz. Data were recorded in two 50 MHz 
wide channels centered at 1.385 and 1.465 GHz, and 
4.835 and 4.885 GHz, respectively. In the 1.4 GHz image 
(2020 s exposure; off-source rms noise = 0.043 mjy/bm) 
shown in Figure [1] we discovered a 3.6" distant knot 
at a position angle (PA) of —66° with a flux density 
of 1.4 mJy (10% uncertainties are assumed in the ra- 
dio measurements for the knot). Because no counterpart 
was detected in the shallower 4.9 GHz image (710 s), we 
obtained deeper (1 hr exposure; rms = 0.022 mjy/bm) 
follow-up B-array observations at this frequency on 2008 
Jan 13 (program S8723) matching the resolution of the 
A-array 1.4 GHz discovery data. The knot is confirmed 
at 4.9 GHz with a flux density of 0.41 mJy and the resul- 
tant 1.4 — 4.9 GHz radio spectral index is a r = 1.0 ± 0.1. 
Model-fitting the radio nucle us and knot in the 1.4 GHz 
(u,v) data using DIFMAP (jShepherd et alj|1994h , we 
found that both features are unresolved at ~ 1" reso- 
lution and set this as the upper limit to the size of the 
knot. 

Our Chandra observation of GB 1428+4217 was ob- 
tained on 2007 Mar 26 (obsid 7874) as part of a small 
snapshot program targeting four z > 3.6 flat-spectrum 
radio quasars with detected radio jets in the aforemen- 
tioned VLA survey (jCheung et alj|2005h . Of these four 
targets, the X-ray jet detections of GB 1428+4217 and 
one other target (PKS 1418-064 at z = 3.7), were ini- 
tially reported in iCheung et all (|2008t ). The full results 
of the survey will be reported elsewhere together with the 
Chandra/YLA results for an additional seven z = 2 — 3 
quasars with radio jets obtained as part of other pro- 
grams. 

In the 11.7 ks Chandra observation of GB 1428+4217, 
we used the nominal aim-point of the ACIS-S3 chip and a 
l/8th sub-array mode (0.4 s frame time) in order to mit- 
igate pileup of the nucleus. We reprocessed the data by 
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Fig. 2. — Chandra ACIS-S 0.5 — 7 keV spectrum of the quasar 
nucleus in the (top) panel. The absorbed power-law model overlaid 
in red, and the deviations from the model in the (bottom) panel. 
The data were grouped in bins with a minimum of 10 counts per 
bin. 

using chandrajrepro script in CIAO 4.4 (|Fruscione et al.l 
2006) and assigned the most recent instrument calibra- 
tion available in CALDB 4.5. The script also runs 
acis_process_events which applies the sub-pixel algo- 
rithm and provides the data with the best angular res- 
olution required by our analysis. With the correction 
for 9.4% deadtime, we have 10.6 ks effective exposure 
time on the source. The astrometry was set by adjust- 
ing the X-ray core posi t ion to that of the VLBI radio 
position from iFev et all (pOOl . R.A. = 14 h 30 m 23 s .742, 
Decl. = +42°04'36.49" (J2000.0). 

More than 2200 counts were detected from the quasar 
X-ray core, allowing for spectral analysis. Model fitting 
of the quasar data was performed with a 1.5" radius cir- 
cular aperture. A pie region excluding the jet was used 
to determine the background. The spectrum was well 
fit with an absorbed power law model (Figure [2]) , with 
fixed Galactic absorption, Nr. The Chandra derived pa- 
ra meters are cons i stent with the XMM results published 
bv lWorslev et alj (|2006l) . with the power-law spectral in- 
dex slightly smaller in our Chandra spectrum, which may 
be due in pa rt to pileup. Applying the jdpileup model 
(jDavisI [200 If ), the core is only 2.7% piled up and the 
resulting parameters arc essentially unchanged, but the 
statistical uncertainties increased (see Table [1]). 

X-ray emission coincident with the 3.6" distant radio 
jet knot is apparent in the smoothed Chandra-ACIS im- 
age (Figure [T]) . The charge transfer readout streak is at 
PA = 43°, i.e., almost perpendicular to the jet axis, so 
does not contaminate the knot's emission. For the knot, 
we defined an elliptical aperture region giving 20.3 ± 4.8 
net counts in the 0.3 — 7 keV range. The spectrum of this 
emission extends up to a rest frame energy of ~ 6.5 (I+2) 
keV = 37 keV, and is quite soft with 14.9 net counts 
between 0.5 — 2 keV and only 4.4 between 2 — 7 keV. 
The fluxes resulting from a spectral fit assuming a sin- 
gle power law with a x = 0.7 frozen and Galactic Nn 
are presented in Table [TJ For the given flux of the jet, 
its isotropic luminosity as measured in the observer rest 
frame in the 0.5—10 keV band would be 4. Ox 10 45 ergs -1 . 
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TABLE 1 
GB 1428+4217 X-ray Spectral Fits 



Component Model 


«x 


N H ( 2 =4.72) 


F(0.5-2 keV) 


F(2-10 keV) 


(1) 


(2) 


(3) 


(4) 


(5) 


Core Power-Law 

Core Power-Law with pileup 


04 o+0.25 


< 3.6 

< 4.9 


5.5 ±0.5 
6.1 ± 0.7 


16.7±S;| 


Jet Power-Law 


0.7 (fixed) 




5.7±H x 10- 2 





(1) Component and source model. All fits assume Galactic Nh = 1.4 xlO 20 cm~ 2 fixed (Dickcv & Lockman 1990). Uncertainties are 90% 
for one significant parameter and upper limits are quoted at 3tr. 

(2) We use the definition of spectral index, a, as F v oc v~ a . 

(3) Intrinsic absorption at 2=4.72 in units of 10 22 cm -2 . 

(4, 5) Unabsorbed flux in the observed energy range in units of 10 -13 erg s~ 1 cm~ 2 . The absorbed F(0.5— 2 keV) values are 14 — 15% 
smaller assuming that the absorber has the Nh value set at the 3tr limit. 



To help constrain the overall spectral energy distribu- 
tion (SED) of the 3.6" knot, we analyzed an archival 
HST WFPC2 image from 1999 Feb 8 (total exposure: 
9.4 ks) with the F814W filter (program 7266). Within 
a circular aperture (r = 0.5") centered on the knot, we 
found no significant excess optical emission with respect 
to the background measured from six adjacent regions. 
From the aperture count rates, we derive an aperture 
corrected fH oltzman et al.iri995| ) upper limit of < 0.13 
/iJy (3a) at 3.74 x 10 14 Hz from the knot. 

Finally, to search for radio jet emission on smaller 
scales, we analyzed VLBI 2.3 and 8.6 GHz data for 
GB 1428+4217 obtaine d in 1998 Feb 1 as part of the 
USNO RRFIE0 project (IFev et al .112001 . The calibrated 
(u, v) data were reimaged and self-calibrated, revealing a 
one-sided jet extending west of the core (Figure [3]). The 
jet is visible in the 2.3 GHz map out to 30 mas (~200 
pc, projected). This jet can be modeled as three knots 
at PA = -96° to -117° with a total flux density of 9.5 
mJy. Utilizing the closest (and brightest) jet knot, we 
measure a jet to counterjet flux ratio defined as the peak 
/ (3x rms noise) = 3.3 mJy / 3 x 0.16 mJy sx 7. Most 
of the pc-scale jet is resolved out in the higher resolution 
8.6 GHz image, where the je t extends ou t to on ly 3 mas 
(2.4 mJy). On smaller scale. IVeres et al.l (|2010D found a 
faint ~ 1 mas extension in VLBA 15 GHz data also in 
the same direction. The pc-scale jet deviates by about 
30° from the large-scale Chandra /VL A knot; however, its 
existence makes the jet interpretation for the kpc-scale 
emission more likely. 

3. DISCUSSION 

The kpc-scale radio and X-ray jet reported here in the 
z = 4.72 blazar GB 1428+4217 is the highest-redshift 
example thus far. Together with the blazar nature of 
its core emission (Section 1), the one sidedness of the 
kpc-scale and ~200 pc scale jet seen in our VLBI images 
imply that the jet is relativistic and probably aligned at a 
small angle, 6, to our line of sight. For 9 20°, the 3.6" 
knot distance from the core corresponds to 24/ (sin 6) 
kpc 70 kpc, deprojected for the detected jet. This 
is far enough from its parent host galaxy that the most 
significant source of seed photons for inverse Compton 
scattering is the cosmic microwave background (CMB). 
This is especially relevant at z=4.72, where the CMB en- 

9 This research has made use of the United States Naval Obser- 
vatory (USNO) Radio Reference Frame Image Database (RRFID). 



ergy density for an observer at rest at the source redshift 
is u* CMB = 4.2 x 10" 13 (1 + z) 4 = 4.5 x lO" 10 erg cm" 3 , 
i.e., 1070 times greater than it is locally. 

The radio and X-ray morphology of GB 1428+4217 
at z = 4.72 is nearly identical to that of the 
z = 4.3 quasar GB 1508+5714 with sim ilar X-ray 
( Siemiginowska et al.ll2003t lYuan et al.ll2003T ) and radio 
(jCheungl 120041: ICheung et al.l |2005[) luminosities in the 
jet. Both show single radio features separated from 
the bright nucleus and display steep spectra (with spec- 
tral indices, a r = 1.0 ± 0.1 and 1.4 ± 0.2, respec- 
tively). Because of poor statistics, the X-ray spectrum 
in GB 1428+4217 is unconstrained and can not be com- 
pared with that of GB 1508+5714 where a x » 0.9 ± 0.4 
was determined. In GB 1428+4217, we calculate the 
ratio / x //r = (^xix) / (fr-Fr) = 205, using monochro- 
matic flux densities at 1 keV and 1.4 GHz. This is larger 
than the z — 4.3 case where f x /fi = 1 58, and larger 
than found in other lower- z examples (cf.. lCheung|l2004l : 
iMassaro et ail 1201 ID . Equivalently, / x // r = 205 corre- 
sponds to a radio to X-ray spectral index, a rx = 0.72. 
Although a rx is just consistent with the limit on the ra- 
dio to optical spectral index, a vo > 0.73, it is smaller 
than the radio spectral index, a r — 1.0, indicating that 
the X-rays are not a simple extension of the radio syn- 
chrotron spectrum. 

The X-ray emission could correspond to either in- 
verse Compton emission of the low-energy electrons in- 
volving CMB target photon field ('IC/ CMB' model; 
iTavecchio et~aT] l2000t iCelotti et all [2001), or an addi- 
tional synchrotron com ponent due to a higher-energy 
electr on population (e.g. JStawarz et a l. 2004; Hardcastlc 
120061 ). Overall, the broad-band spectral indices appear 
consistent with the (1+z) 4 amplification of the CMB en- 
ergy density, / x // r ~ (5/T) 2 u' CMB /u' B oc (1 + z) 4 (6/B) 2 
(with the Doppler factor, 6, the bulk Lorentz factor, T, 
and the CMB and magnetic field energy densities in the 
jet rest frame, u' CMB — T 2 u^, MB and u' B — B 2 /8n 1 re- 
spect ively), as would be expected in the IC/CM B model 
e.g.. jSchwartzl 120021: ICheung 2004). Following [Chgungl 
2001 . we apply the IC/CMB model assuming also 
equipartition with a relativistic electron spectrum with 
power-law slope, p = 2a + 1 = 3, extending down to a 
minimum 7 = 10, in a sphere with an appar ent size of 
1.8 x 1" for the GB 1428+4217 jet knot (see iMarscherl 
1983). For the case where the jet Lorentz factor is set 
equal to the Doppler factor, we derive V = S — 4.7, and 
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Clean ER map. Array: BFGHKLN0P5G BFGKKLMMNOPSG 
J1430+4204 at 2.292 GHz 1998 Feb 10 



Clean RE map. Array: BFGHKLMNOPSG BFGKKLMMNOPSG 
J1430 + 4204 at 8.646 GHz 1998 Feb 10 
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Fig. 3.— VLBI 2.3 GHz (left) and 8.6 GHz (right) images of GB 1428+4217 from reprocessing USNO data HFev et al.ll2"004T) . Note the 
difference in scales. The contours begin at 0.5 (3.87 mas X 3.08 mas) ) and 0.4 mjy/bm (1.09 mas X 0.73 mas), respectively, and increase 
by factors of \/2. The bright core with the one-sided pc-scale jet to the west is apparent. 



a magnetic field in the jet rest frame, B = 35 /iG. These 
parameters are fairly insensitive to the extrapolation of 
the p — 3 slope down to low energies and our conser- 
vative upper limit to the radio knot size. Assuming in- 
stead a typical observed radi o jet spectral index of 0.7 5 
at lower-redshifts (following iKataoka fc Staward 12005). 
and the knot diameter is 0.6 x smaller than the assumed 
1" limit, we obtain B = 21 [iG and S — 6.3. These esti- 
mates assume a single X-ray /radio emitting zone applies 
to our ~ 1" (6.6 kpc) angular resolution element. 

Similarly applying the IC/CMB model to t he other 
high-redshift X-ray jet, GB 15 08+4714 (z = 4.3:lCheund 
[200l and 1745+624 (z = 3.9: ICheung et al.ll2006D . thev 
derived 5 ~ 3 — 5. These values are smaller than the 
typical o nes (<5 ~ 5 — 10) determin ed for lower-redshift 
jets (e.g.. IKataoka fc St awarz 2005). Although the num- 
ber of quasar X-ray jets found at such very high-redshifts 
(z ~ 4—5) is still small, the slower large-scale jets implied 
in our analysis for the current sample seem to indicate 
that they are either: (1) intrinsically less relativistic, or 
(2) decelerate more rapidly out to ~ 10's — 100's kpc- 
scales than their lower-redshift counterparts. In the first 
case, this would be consistent with the recent finding that 
the relative abundance of high-power blazars relative to 
the parent population of radio-loud quasars and radio 
galaxies decreased substantially above z ~ 3, likely as a 
result of a decrease of the average bulk Lorentz factor 
of bl azar jets in the early Universe (see IVolontcr i et al.l 
2011). The latter possibility could be supported by the 
fact that the environments of high-redshift (z > 2) radio 
sources are believed to reside in more inhomogeneous and 
multi-ph ase environments than in n earby radio sources 
(see e.g.. lReeslH98llDe Youm3 l2006). and this may man- 
ifest in the slower large-scale jets due to an enhanced 
entrainment of the ambient gas when the jet propagates 
through the host galaxy. 

An argument put forth against the IC/CMB model 
comes from observations of the unresolved (at Chandra 



resolution) quasar "cores" at high-redshift. In this inter- 
pretation for the kpc-scale emission, one may expect that 
the X-ray emission in these cores could also be enhanced, 
as they contain an unresolved portion of the jet (which 
we know are highly relativistic from, e.g., VLBI superlu- 
minal motion studies). This could manifest in different 
X-ray core properties at high-z (X-ray, and optical-to-X- 
ray spectra), but Chandra studies thus far sh ow no such 
trends with redshift (e.g.. lLopez et al.ll2006[ ). However, 
these analyses are hampered by the fact that most of 
the sources analyzed so far in this context are only mod- 
erately radio-loud, and in general, only a small fraction 
of jets in radio-loud quasars extend to k pc-scales (e.g., 
IBridle & Perk^fll^lLiu" fc Zhandl2002T) . 

The radio structures of higher-red shift samples a re now 
being investigated in more detail (|Gobeilld)20ill Gob- 
eille et al., in preparation), and will pave the way for 
further Chandra X-ray studies. X-ray observations of a 
larger sample of kpc-scale radio jets over a broad redshift 
range (z > 2) should provide further elucidation of the 
correct X-ray emission mechanism which will give us in- 
sight on the deeper physics contained in the data. Such 
high redshift systems allow us to study jets which result 
from the "earliest" actively accreting black hole systems. 
As most detections are currently at z £ 2, we have begun 
to study more distant examples and these results will be 
presented in a future paper. 
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